Eur J Appl Physiol (2003) 88: 317-338 
DOI 10.1007/s00421-002-0676-3 


REVIEW ARTICLE 


Alois Mader 

Glycolysis and oxidative phosphorylation as a function of cytosolic 
phosphorylation state and power output of the muscle cell 


Accepted: 12 June 2002 / Published online: 22 November 2002 
© Springer-Verlag 2002 


Abstract A mathematical description of the regulation 
of ATP production in muscle cells is presented whereby 
the activity of OxP can be calculated as a function of (1) 
free [ADP] as the substrate and (2) a second driving 
force ®AG (kilojoules per mole) resulting from the 
difference of free energy AG ox , ap (kilojoules per mole) - 
AGATP,cyt (kilojoules per mole). In turn, the term AG oxap 
results from the proton motive force and the generation 
of ATP in the matrix space including the ATP-ADP 
exchange, whereas the phosphorylation state of the 
CHEP-sytem is described by AGATP.cyt- Regulation of 
glycolysis is calculated as a function of free [ADP] and 
[AMP] at the level of PFK. The PFK is inhibited by a 
decreasing pH resulting from lactate accumulation. The 
ATP/PCr equilibrium of the CHEP-system is calculated 
by algebraic equations. The dynamic behaviour of the 
metabolic control of ATP production as a function of 
ATP consumption is calculated by a system of two 
lst-order non-linear differential equations, including a 
time delay considering oxygen transport. Lactate 
distribution and elimination is calculated using a two- 
compartment model with an active lactate producing, 
and a passive, space including lactate elimination by 
combustion. The simulation of the dynamics of energy 
metabolism of muscle cells is performed by the stepwise 
solution of the differential equations with a 5th-order 
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Runge-Kutta-Fehlberg-routine. Examples of various 
applications of the simulation of the dynamics of energy 
supply demonstrate the qualitative and quantitative 
congruence to the behaviour of metabolic processes in 
experiments during rest, exercise and recovery. 

Keywords Regulation of oxidative phosphorylation 
and glycolysis • Mathematical model • Computer 
simulation of energy metabolisms • Cytosolic ATP/PCr 
equilibrium 


Symbols and abbreviations 


^^ATP,cyt 

Current free energy of the CHEP sys¬ 
tem (at rest -58 and -65 kJ-mol ') 

^^ATP,mit 

Current free energy of the inner mi¬ 
tochondrial ATP/ADP system build up 
by Ap and the FoFi-ATPase 

A^ox,ap 

Current free energy build up by AGATP.mit 
and the ATP 4 /ADP~ 3 translocase 

^^ox,const 

-72.5 kJ-mol ’, Arbitrarily set maxi¬ 
mum of free energy of AG ox>ap used in 
the simulation model 

o 

< 

Standard free energy of the [ATP]/[PCr] 
system, -30.6 kJ-mol 1 

Ap 

Mitochondrial proton motive force 
(pmf), about -220 mV, equivalent to 
about -21.7 kJ-mol 1 generated by the 
respiratory chain 

Ap' 

pmf expressed in kilojoules per mole 

‘LAG’adp 

Force equivalent of free [ADP] result¬ 
ing from AG ox>ap -AG A TP,c y t, being the 
driving force necessary to sustain the 
flow of ATP from inside the mito- 
chondrian to the cytosol 

Ai ]/ 

Mitochondrial transmembrane electri¬ 
cal potential, about -160 mV 

Z[A] 

Sum of concentrations of cytosolic ad¬ 
enosine phosphate compounds (S[A] = 
[ATP] + [ADP] + [AMP]), millimoles 
per kilogram 
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t F02 

used as a refence for an approximate 
characterisation of the response delay 
of a non lst-order system by the time 
constant (or delay time) of a first order 
system. The time behaviour of a lst- 
order system is described by: a = e 0 -( 1- 
e~ //T ) where a is the response amplitude 
at time t after the start, and e 0 is the 
final amplitude approached after an 
infinite time t. The r = e 0 ( 11 /e) (ap¬ 
proximately 0.666-e o ) 

Z[C] 

Sum of concentrations of compounds 
available to store phosphate bond en¬ 
ergy in the cytosol (S[C] = [PCr] + [P] ; ), 
millimoles per kilogram 

[ADP] 

Concentration of cytosolic adenosine 
diphosphate, millimoles per kilogram 
wet mass, or millimoles per litre 

[ADP] mit 

Concentration of inner mitochondrial 
adenosine diphosphate, matrix space, 
millimoles per litre 

[AMP] 

Concentration of adenosine mono¬ 
phosphate, millimoles per kilogram wet 
mass, or millimoles per litre 

[ATP] 

Concentration of cytosolic adenosine 
triphosphate, millimoles per kilogram 
wet mass, or millimoles per litre 

[ATP] mit 

ATP 4 /ADP 3 - 

Concentration of the inner mitochond¬ 
rial ATP, matrix space, millimoles per 
litre 

translocase 

protein complex located in the in¬ 
ner mitochondrial membrane which 
transfers ATP (built up by the F 0 F r 
ATPase) from matrix to the cytosol in 
exchange for ADP. The ATP" 4 and 
ADP" 3 are transferred in ionic form; so 
ATP is 4 times and ADP is 3 times 
negatively charged 

bpow 

Rate of GP consumption per watt, 
millimoles per kilogram per second per 
watt 

b v ,i a 

Transformation coefficient of v [a into 
v GP vi a , about 1.4 mmol lactate-mmol 1 
GP 

t>F02 

Transformation coefficient of F() 2 into 
v G p,F02, a , a bout 4.3 ml0 2 -mmol 1 GP 

Ca 2+ 

Calcium ion 

CHEP 

Cytosolic high energy phosphate- 

CK 

Creatine kinase 

C P /•; 

Arbitrary conversion factor 

[Cr] 

Concentration of creatine, millimoles 
per kilogram muscle 

CS 

Computer simulation 

d bu ff 

Coefficient equal to the mean buffering 
capacity of skeletal muscle 


d[GP]/d? Rate of change of the sum of 

[ATP] + [PCr], millimoles per kilogram 
per second 

d[GP]j>o 2 , a /df Rate of change of instantaneous oxi¬ 
dative GP production 

d[la ] b /dt Rate of change of lactic acid concen¬ 

tration (millimoles per litre per second) 
in the passive lactate producing com¬ 
partment (blood) 

d[la ] m /d/ Rate of change of lactic acid concentra¬ 

tion (millimoles per litre per second) in 
the active lactate producing compart¬ 
ment (muscle) 

dv G p,F02,a/dt Rate of change of v G p,no 2 ,a 

dV0 2 Jdt Rate of change of F0 2 , a 

E Energy flow, watts 

FoF^ATPase protein complex located in the inner mi¬ 

tochondrial membrane which harvests 
the energy of 3 protons to build up one 
molecule of ATP. The energy transferred 
by the move of the protons through the 
F 0 proton channel results from the pmf 
build up across the inner mitochondrial 
membrane. The ADP binding of the F r 
complex leads from the open to the low 
state where ADP is weakly bound. The 
move of H + from the cytosolic to the 
matrix side builds up the tense state until 
GATP.mit is achieved and ATP is released 
from the F,-complex 

F ag Factor which adjusts current power 

output [F AG -E(watts per kilogram)] in 
relation to AG atp cyt 

[GP] Sum of concentrations of high energy 

phosphate compounds ([GP] = [ATP] + 
[PCr]), millimoles per kilogram 
[FI + ] Hydrogen ion concentration 

K t Coefficient relating rate of diffusion of 

lactate from muscle in terms of [la 
given in Eq. 31 

k ebox 50% Activity rate constant of oxidative 

lactate elimination by combustion, as a 
function of lactate concentration 
Km(0.5) Common expression for the 50% rate 

constant of OxP. Equivalent to v ss / 
Vax 0.5 

ksj 50% Activation rate constant of Eq. 13 

of OxP related to F0 2max 

ks 1DG From -0.65 x 10" 4 to -0.95 x 10 4 -moF 1 

free [ADP] per joules per mole free 
energy of the CHEP system 
ks'j Hill coefficient of 50% activity of Eq. 9 

of OxP related to F0 2max 

ks 2 50% Activation rate constant of glyc¬ 

olysis mainly due to PFK activation, 
[ADP] and [AMP] as activating factors 
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kS3 

50% Inhibition rate constant of glyc¬ 

R 


olysis due to PFK inhibition by de¬ 

RQ 


crease of pH m 

T 

knc>2 

reciprocal of ty 02 , per second 


[la] 

Concentration of lactate, millimoles per 
litre 

r V02 

[la ] b 

Concentration of lactate in blood, mil¬ 



limoles per litre 

v GP,F02,a 

[la ] m 

Concentration of lactate in muscle, 
millimoles per litre 


[la ]ss 

Steady state concentration of lactate 
(millimoles per litre) during an equilib¬ 
rium of lactate production and oxida¬ 

v GP,v,la 


tive lactate elimination (0 = Vi a>ss ^ pH - 

Vla.ox) 

V GP,£ 

Mj 

Hydrogen ion dependent creatine kin¬ 



ase equilibrium constant 

v GP,rest 

m 2 

Hydrogen ion independent creatine 
kinase equilibrium constant 



1.66 x 10 9 mol _1 

Via 

m 3 

Adenylate kinase equilibrium constant, 


Mg 2 + 

range 0.85-1.15-mof 1 

Vla.max 

Magnesium ion 


mvalF 1 

Expression for bases or acids of sub¬ 



stances having an ionic dissociation in 
H 2 0 whose concentration is given in 
mmoH . The unit mval expresses the 
ionic strength of the acid base status of 

V|a.max,pl 1 


blood and tissues 

Vla.ox 

NAD + 

Nicotinic adenine dinucleotide (oxi¬ 
dized form) 

NADH 

Nicotinic adenine dinucleotide (reduced 
form) 

Vla.ox.b 

nox 

Exponent of the Hill equation for cal¬ 



culating the rate of OxP as a function of 
free [ADP] 

Vla.ox.m 

OxP 

Oxidative phosphorylation 


pco 2 

Partial pressure of carbon dioxide (mil¬ 
limetres of mercury) at the tissue level 

Vla.ss 

[PCr] 

Concentration of phosphocreatine 



(PCR), millimoles per kilogram muscle 

Vla.ss.pH 

PDH 

Pyruvate dehydrogenase 


PFK 

Phosphofructokinase 


pH 

~log 10 [H + ] 

C mit 

pH m 

Intracellular cytosolic pH 

I/e 1 

[P]i 

Concentration of free inorganic phos¬ 
phate. millimoles per kilogram muscle, 
or millimoles per litre 


PK 

Pyruvate kinase 


pmf 

Proton motive force 


31 P-NMR 

3 'Phosphorus nuclear magnetic reso¬ 
nance spectroscopy 


P/O quotient 

Relationship between the gain of high 



energy phosphate bonds and the H 2 0 
produced. Typical value about 2.6 

1.35 


The gas constant, 2.0 x 10 3 kcal-moF 1 
Respiratory quotient 
The absolute temperature Kelvin 
(C + 273) 

Time constant(seconds)of the rise of 
V0 2 from the beginning of exercise, 
and during a period of constant load 
Rate of oxidative rephosphorylation 
expressed in terms of [GP], millimoles 
per second per kilogram 
Rate of glycolytic rephosphorylation 
expressed in terms of [GP], millimoles 
per second per kilogram 
Power output expressed in terms of rate 
of GP consumption during contraction, 
millimoles per second per kilogram 
Power output expressed in terms of rate 
of GP consumption at rest, millimoles 
per second per kilogram 
Rate of lactic acid formation, millimo¬ 
les per second per kilogram 
Maximal rate of glycolysis expressed as 
rate of lactic acid formation, millimoles 
per second per kilogram 
pH Dependent maximal rate of glycol¬ 
ysis, expressed as rate of formation of 
lactic acid, millimoles per second per 
kilogram 

Rate of lactate elimination by oxida¬ 
tion, millimoles per minute per kilo¬ 
gram 

Rate of lactate elimination by oxidation 
(millimoles per minute per kilogram) in 
the passive compartment (blood) 

Rate of lactate elimination by oxidation 
(millimoles per minute per kilogram) in 
the active compartment (muscle) 
Steady-state rate of gross lactic acid 
(pyruvate) formation at pH above 7.4, 
mmol-s _1 -kg' 1 

The pH dependent steady-state gross 
rate of formation of lactic acid, milli¬ 
moles per second per kilogram 
Mitochondria volume, millilitres 
Quotient %V m /%V b, the volume of 
active lactate producing space (m, 
muscle, about 30% of total body vol¬ 
ume) divided by the volume of the 
passive space (b, blood, about 15% of 
total body volume) available for lactate 
distribution during work and rest. Ac¬ 
tive and passive spaces together con¬ 
stitute between 35% and 50% of total 
body volume 

Factor which relates muscle water space 
to active muscle mass (about 0.72) 
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vo 2 

Oxygen uptake, millilitres per minute 
or, millilitres per minute per kilogram 
body mass 

vo 2 , d 

Current V0 2 , millilitres per second per 
kilogram muscle 

I 1 f\a.i V sl 

Current V0 2 at rest, millilitres per sec¬ 
ond per kilogram muscle 

vo 2max 

Maximal oxygen uptake, millilitres per 
second per kilogram muscle or body 


mass 

V 02max,ap 

Simulated exponential increase in oxy¬ 
gen uptake, approaching the actual 
oxygen uptake, millilitres per second 
per kilogram muscle or body mass 

vo 2 ' SS 

Steady state of V0 2 , millilitres per sec¬ 
ond per kilogram muscle mass wet 
weight 

watts 

Power, joules per second 


Introduction and the general aspects of the regulation 
of the ATP producing reactions 

The maintenance of a state of high phosphorylation of 
the CHEP-system and maintenance of the metabolic 
rate of ATP production is essential for the life of a cell; 
this is achieved using feedback control which adjusts 
the activity of the rephosphorylating processes to the 
instantaneous rate of dephosphorylation. There has, 
however, been disagreement concerning the explanation 
of the regulation of OxP in vivo for about 20 years. 
The regulation of the activity of OxP in vivo at a cel¬ 
lular level is believed to be different from that occur¬ 
ring in vitro (Balaban 1990; Brown 1992; Hochachka 
and Matheson 1992; Hochachka and McLelland 1997). 
The in vivo experimental results show that there is a 
gap between the flow rate of OxP, as obtained by the 
Michaelis-Menten (MM)-type kinetics as a function of 
the calculated free [ADP], and the experimental ob¬ 
servations of the rate of OxP (Balaban and Heineman 
1989; Connett 1988; Connett and Honig 1989; Hein¬ 
eman and Balaban 1990; Kantor et al. 1986; Katz et al. 
1988a, 1989; McMillin and Pauly 1988). Many authors 
therefore have tended to neglect the role of free [ADP] 
as one of the main regulating factors, in favour of 
other activation mechanisms (Balaban 1990; Balaban 
and Heineman 1989; Brown 1992; Conley et al. 1997; 
Connett 1988; Connett and Honig 1989; Heineman and 
Balaban 1990; Hochachka and Matheson 1992; Ho¬ 
chachka and McLelland 1997; Kantor et al. 1986; Katz 
et al. 1989; McMillin and Pauly 1988). These other 
potential regulatory mechanism have generally been 
believed to depend on Ca 2+ activation of NAD + de¬ 
pendent dehydrogenases (Territo et al. 2000) and of 
other mitochondrial key enzymes such as PDH 
(Hansford and Zorov 1998). The rise of Ca 2+ caused 
by muscle contraction has therefore been assumed to 
enhance the up-regulation of OxP (Territo et al. 2000). 


Finally, some authors have proposed the inclusion of 
known potentials - the AGATP.cyt and the Ap - in the 
regulation scheme (Barstow et al. 1994; Fontaine et al. 
1997). 

This paper is restricted to a basic quantitative 
analysis of the relationship of the cytosolic phospho¬ 
rylation state to the activity of OxP and glycolysis. It is 
the aim of this paper to demonstrate that the main 
control of the rephosphorylating processes in the cells 
can take place at this level, so that all contradictions 
discussed above regarding the regulation mechanism of 
OxP can be resolved. It is also intended to demonstrate 
that mathematical models derived from the above 
analysis can be used as tools to simulate the complex 
dynamic behaviour of metabolic reactions such as OxP 
and glycolysis, thus contributing to the development of 
a detailed, consistent, theory of the functions involved. 


Basic aspects of a mathematical description 
of the feedback control of OxP and glycolysis 
by the phosphorylation state of the CHEP system 

The feedback control of glycolysis and OxP requires that 
the phosphorylation state of the CHEP system is tightly 
coupled to the rate of glycolysis and OxP, so that a 
certain degree of the phosphorylation deficit activates 
OxP and glycolysis, to compensate for an increased 
rate of dephosphorylation, at rest or during exercise. 
Feedback control systems of this type are dynamic sys¬ 
tems which exhibit two stages: 

1. A steady state which is achieved when the rate of 
ATP consumption is equal to the rate of ATP pro¬ 
duction. At a metabolic steady state the CHEP sys¬ 
tem is at equilibrium; the metabolic rates and the 
main parameters of the internal environment are 
constant. 

2. Transient states caused by sudden changes of the 
ATP consumption rates. This may lead either to a 
new steady state, or to a stage of a non-steady state 
which can be maintained only for a limited time and 
within certain limits of the deviation of the inner 
physico-chemical environment. 

Feedback control systems describing the dynamics of 
the metabolic reactions at the cellular level cannot be 
established and described without the analysis of the 
time behaviour of the main parameters of the system 
(Funk et al. 1990; Mader and Heck 1996; Mader 1998; 
Meyer 1988). The mathematical description must also 
include a sufficient set of parameters and equations 
which determines its behaviour. Then the translation to 
a computer model allows the comparison of the exper¬ 
imental data generated by the simulation of the system's 
dynamic to the measured parameter of the real system as 
obtained by a given experimental procedure. The 
mathematical description of the elements of the meta¬ 
bolic control systems must be based on a sufficient set of 
equations governing the parameters of the system and 
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determining its behaviour. This includes the calculation 
of: 

1. The equilibrium constants of the several reactions 
among the components of the CHEP system, in¬ 
cluding the calculation of AGATP.cyt (kilojoules per 
mole). 

2. The steady state activity characteristics of glycolysis 
and of OxP. These will be described as functions of 
the phosphorylation state of the CHEP system, the 
known activators and the equivalent of Ap' (kilojo¬ 
ules per mole). This will be expressed in terms of 
AG ox>ap (kilojoules per mole) which results from the 
inner mitochondrial ATP formation by the F 0 Fi- 
ATPase complex and the ATP-ADP-transfer by the 
ATP 4 /ADP 3 -translocase. An excellent explanation 
of these functions is given in Ferguson (2000) 

3. The influence of some main parameters of the inner 
chemical environment of the cell on the chemical 
equilibrium of the CHEP system and on the activity 
of OxP and glycolysis. 

On these bases, a consistent theoretical description of 
the main regulation of OxP and glycolysis will be de¬ 
veloped, leading to the development of mathematical 
models of differing accuracies. These will be shown to 
be essential for an extended quantitative recalculation of 
the steady state and of the dynamic behaviour of the 
main parameters of this system, and to be related to all 
possible experimental situations. This will be feasible 
only by using differential equations for the calculation of 
the parameter of the system as a function of time and 
“functional load” or the disturbances of the system 
(Mader and Heck 1986, 1996). 

Calculation of the equilibrium between PCr 
and adenosine phosphates 

The equilibrium of the cytosolic high energy phosphate 
system results from the CK reaction: 

A. [ATP] + [Cr] <=====> [ADP] + [PCr] + [H + ] 
including the adenylate kinase reaction: 

B. 2[ADP] <=====> [ATP] + [AMP] 

However, reaction (A) is not a real reaction. Ac¬ 
cording to Golding et al. (1995) and others 
(McGilvery and Murray 1974) the real reactions 
which are summarized under ( A) consist of a set of 
about 12 reactions including free [P] ; , H + and 
Mg 2+ ions. In contrast to the common interpreta¬ 
tion of reaction (A), it can be assumed that [Cr] 
functions as a template for phosphate binding and 
the CK reaction can be rewritten as a mirror reac¬ 
tion where the free cytosolic [P]; replaces [Cr]: 

C. [ATP] + [P]i < = = = = = > [ADP] + [PCr] + [H + ] 

This takes into account the fact that, under condi¬ 
tions at rest, the content of free [Cr] can vary between 
5 and 15 mmol-kg 1 wet muscle tissue (Arnold et al. 
1984; Conley et at. 1997) but free cytosolic [P] ; is kept 


at a relative low concentration of about 2 to 4 mmo¬ 
l-kg 1 in normal skeletal muscle (Arnold et al. 1984; 
Barstow et al. 1994; Conley et al. 1997; Connett and 
Honig 1989; Kushmerick 1987; Kushmerick et al. 1992; 
Nicholls and Ferguson 1997; Nioka et al. 1992; Taylor 
et al. 1986; Thompson et al. 1995; Williamson 1979; 
Wilson 1994). This applies also to muscles with a high 
oxidative capacity, such as the heart (Balaban 1990; 
Heineman and Balaban 1990; Kantor et al. 1986; Katz 
et al. 1988a, 1989; Lewandowski et al. 1987). This in¬ 
dicates that the ATP-PCr equilibrium in vivo is de¬ 
termined by AGATP.cyt °f ^58 to -65 kJ-mol 1 
calculated by Eq. 4 which rules the [PCr]:[P] ; ratio in¬ 
dependently of the excess of [Cr]. Therefore, the free 
energy of the CHEP system AGATP.cyt is calculated 
using [P] ; rather then [Cr]; it could be rewritten also as 
a function of the ratio [PCr]:[P] ; if the limited amount 
of free cytosolic [P] ; under conditions of rest is taken 
into account (Barstow et al. 1994). Without this pre¬ 
caution, there is a theoretical inconsistency between the 
calculation of free [ADP] and AGATP.cyt on the one 
hand and the free cytosolic [P]i and [Cr] on the other. 
This inconsistency might be one of the causes limiting 
the understanding of the experimentally known facts 
regarding the regulation of OxP. Regardless of the 
excess of [Cr], the free cytosolic [P]j under conditions 
of rest is kept within the range, determined by 
AGATP.cyt, °f about -58 to -65 kJ-moF 1 as calculated 
using Eq. 4 (Arnold et al. 1984; Balaban 1990; Barstow 
et al. 1994; Brown 1992; Heineman and Balaban 1990; 
Kantor et al. 1986; Kushmerick 1987; Kushmerick et al. 
1992; Nicholls and Ferguson 1997; Nioka et al. 1992; 
Taylor et al. 1986; Thompson et al. 1995; Williamson 
1979; Wilson 1994). Since free Mg 2+ concentration 
does not influence the CK-equilibrium much, provided 
it is 1.0 nnnol-1 1 or more (Golding et al. 1995; Veech 
et al. 1979), its influence will not be taken into ac¬ 
count. 

The chemical equilibrium between the components of 
the CHEP-system is established within a few millisec¬ 
onds, so, a change of the [ATP]:[ADP] ratio, brought 
about by dephosphorylation or rephosphorylation pro¬ 
cesses, leads to an immediately re-established ATP-PCR 
equilibrium, as indicated by the saturation transfer 
measurement between [ATP]:[ADP] and [PCr]:[P]; on 
living muscle using 31 P-NMR (LeRumeur et al. 1997). 
Thus, it can be assumed that under all conditions in 
living cells, the CHEP system is near to a simple chem¬ 
ical equilibrium (Arnold et al. 1984; Barstow et al. 1994; 
Kushmerick 1987; McGilvery and Murray 1974; Taylor 
et al. 1986; Thompson et al. 1995; Wieseman and 
Kushmerick 1997). 

The equilibrium between [PCr]:[P]; and the adenosine 
phosphates ([ATP], [ADP] and [AMP]) including the 
influence of [H + ] can be calculated at different levels of 
accuracy (Golding et al. 1995; McGilvery and Murray 
1974). The assumption that reaction (C) is the mirror 
reaction of (A), coupled with reaction (B) provides the 
simplest algebraic solution in the calculation of the 
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ATP-PCr equilibrium during dynamic changes in the 
concentration of the CHEP system (see Eqs. 6, 7 and 8). 
According to reaction (C) [ATP] and [PCr] at equilib¬ 
rium are described by: 

1) [ATP]/[ADP] = M, x [PCr]/[P]j (1) 

with M| equal to the apparent equilibrium constant. 

Since the ATP-PCr equilibrium depends on pEl ra , the 
apparent CK equilibrium constant Mi is given by the 
product M 2 x [H + ] [M 2 = 1.66 x 10 9 -mol 1 (Veech et al. 
1979)]. Rearranging Eq. 1 the level of free [ADP] can be 
calculated as: 

2) [ADP] - [ATP] x [P] ; /(M 2 x [H + ] x [PCr]). (2) 

Free [AMP] can be calculated from [ADP] using the 
following equation (see reaction B): 

3) [ATP] x [AMP] = M 3 x [ADP] 2 (3) 

where M 3 is the equilibrium constant (from approxi¬ 
mately 0.85-1.05) (Krause and Wegener 1996; 
Newsholme and Start 1973). According to Eq. 3 [AMP] 
is related to [ADP] 2 ; therefore the product [AMP] x 
[ADP] is approximately equivalent to [ADP] 3 (Mader 
et al. 1983; Mader and Heck 1986, 1996). 

The free energy of the CHEP system (AGATP.cyt) is 
determined by the quotient of [ATP]/([ADP] x [PJ), 
which represents the state of phosphorylation and can 
be calculated by: 

4) AG A TP,cyt(kJ mol 1 ) =AG 0 - R x T x ln([ATP]/ 

([ADP] x [P] ; )) (4) 

with AG o = -30.6 kJ-mol '. 

If it is assumed that S[A] = [ATP] + [ADP] + [AMP] 
and if S[C] = [PCr] + [P] ; (because free [Cr] is normally 
present in excess over [P] ; ) this system represents the 
amount of substrate available for energy storage. If, as a 
first approach, S[A] and S[C] are assumed constant, the 
phosphorylation deficit of the CHEP system can be ex¬ 
pressed either by the difference S[A]-[ATP] or by the 
difference S[C]-[PCr]. This allows the construction of an 
algebraic solution to calculate free [ADP] as a function 
of the difference S[C]-[PCr], including the influence of 
[H + ] on the [ATP]-[PCr] equilibrium, where S[C] rep¬ 
resents the sum of the active compounds [PCr] + [P] ; 
related to true conditions of rest. 

Furthermore, since the equilibrium assumption is 
also valid for AG A TP,cyn this allows a direct calculation 
of AG A TP.c y t from the quotient [PCr]/(S[C]-[PCr]). Sub¬ 
stitution of [ADP], as given by the right side of Eq. 2, in 
Eq. 4 the range of AGATP.cyt of the CHEP system can be 
obtained from the [PCr]/[P] ; quotient, including the in¬ 
fluence of [H + ], if it is assumed that [P] ; is the chemical 
active factor. Introducing the above definitions ([AT- 
P] + [ADP] + [AMP] = S[A] and [PCr] + [P] ; = S[C]) then 
the phosphorylation state of the CHEP system is indi¬ 
cated by the common term S[C]-[PCr] of the following 
Eqs. 5, 6, 7 and 8 (Mader and Heck 1996): 


5) AG A TP,cyt(kJ ■ mol *) =AGo — R x T 

x ln{M 2 x [H+] x [PCr]/(S[C] - [PCr]) 2 } (5) 

6) [ADP] = S[A]/{{(M 3 x (S[C] - [PCr])/M 2 

x[H+] x [PCr]} 

+{M 2 x [H+] x [PCr]/(S[C] - [PCr])}+l} 

( 6 ) 

7) [ATP] = M 2 x [H+] x [PCr] x [ADP]/(S[C] - [PCr]) 

( 7 ) 

including the calculation of free AMP which is: 

8) [AMP] = S[A] - [ATP] - [ADP], (8) 

The use of Eqs. 5 , 6, 7 and 8 for calculating the equi¬ 
librium of the high energy phosphates in the cytosol has 
the advantages that the phosphorylation deficit (= S[C]- 
[PCr]) and AGATP.cyt of the CHEP system are calculated 
on a unified basis and the relationship between [P]; and 
[PCr] can be determined by 31 P-NMR-spectroscopy with 
the possibility of calculating free cytosolic [ADP] as well 
as AGATP.cyt using the same set of experimentally mea¬ 
surable values. 

The interplay of the different variables of the system 
as predicted by Eqs. 5, 6, 7 and 8 is displayed in Figs. 1 
and 2 as a function of the sum [P] ; + [ADP] + [AMP]. 
The above description of the equilibrium conditions of 
the components of the CHEP system, at the simplest 
possible level, also provides the basis for the calculation 
of the activity of glycolysis and oxidative phosphoryla¬ 
tion as a function of the phosphorylation state of the 
CHEP system. This problem will be tackled in the next 
section. 

Calculation of the steady-state activation of OxP 

If the oxygen supply at the tissue level is sufficient, the rate 
of OxP (and as a consequence the rate of V0 2 by mito¬ 
chondria) is to a large extent determined by free cytosolic 
[ADP] as substrate. According to Chance and Williams 
(1955) the in vitro ratio v vss /v max = 0.5 is achieved at about 
0.035 mmol! 1 free [ADP]. The steady-state activity of 
OxP expressed as rate of V0 2ss P er kilogram muscle may 
be calculated using the Hill equation (mathematical 
description of an MM-kinetic having an exponent greater 
than 1: for details see Newsholme and Start 1973) as a 
function of free [ADP] and V0 2max (Mader et al. 1983; 
Mader and Heck 1986. 1996) from: 

9) C0 2jSS =E0 2 , max /(l + ks' 1 /[ADP] nox ) (9) 

The coefficient ks'i is the 50% activation rate constant 
related to the exponent nox. This exponent must be 
greater than 1.0, since using the MM-kinetic with an 
exponent nox =1.0 it is mathematically impossible to 
calculate an appropriate activation of V0 2 as an exclusive 
function of [ADP]. A reasonable value for the exponent 
nox can be assumed to be in the range of 1.4 or higher, to 
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Fig. 1 Phosphorylation relationships according to McGilvery and 
Murray (1974), on a linear scale. The free energy AGatp,cj 7 
(kj-mol -1 ) is calculated using Eqs. 4 or 5, for S[A] = 7.0 mmol-kg" , 
S[C] = 28 mmol-kg 1 muscle, CK equilibrium constant 
Mj = 1.66T0 9 mor' and M 3 =1.05-mor 1 (McMillin and Pauly 
1988; Taylor et al. 1986), for intracellular pH values from 6.2 to 
7.4 in steps of 0.2. The state of phosphorylation for a well- 
oxygenated muscle at rest (pH = 7.0) is in the range indicated by 
Rest. Fatigue is related to pH = 6.4 and [PCr] close to 0 and is 
indicated by Fatigue. In conditions of rest the free energy AG ATP cyt 
of the CHEP-system is close to AG OXjap (-60 to -65 kJ-mol 1 ) 
generated by the apparent Ap. The constant AGox. const = 
-72.5 kJ-mol 1 is related to the rate of OxP at rest, caused by the 
proton leak. As indicated by the arrow , with increasing dephosph¬ 
orylation of the CHEP system, the potential difference <j>AG 
between AGox c °nst and AGatpCYI increases linearly. For more 
information see text and list of abbreviations 

establish a (slightly) sigmoidal activation characteristic 
(Mader and Heck 1996), (see Fig. 3). Expressed in terms 
of free [ADP] ks'j is assumed to be in the range of 
0.035 mmoll 1 (Chance and Williams 1955; Nioka et al. 
1992). 

The values (ks'x = 0.035 mmol l 1 and nox=1.4) will 
be used, even if recent experiments favour a Hill equa¬ 
tion (Newsholme and Start 1973) with Km(0.5) = 
0.048 mmol l 1 free ADP and an exponent close to 2.0 
(Cieslar and Dobson 2000). Indeed, many experimental 
results are consistent with a 50% activity constants 
Km(0.5) related to free [ADP] in the range between 
0.020 and 0.050 mmol l 1 (Barstow et al. 1994; Kushm- 
erick et al. 1992; Nioka et al. 1992; Thompson et al. 
1995). It should also be noted that V0 2max is related to 
the mitochondria volume of the muscle cell and the P/O- 


quotient in living mitochondria is approximately con¬ 
stant and independent from the rate of OxP; therefore 
the rate of OxP is a linear function of V0 2 , as pointed 
out in detail later (Nicholls and Ferguson 1997). 

The functional origin of the exponent nox >1.0 using 
the half maximal rate constant ks' t approximating 
0.035 mmol l 1 free ADP needs a thorough explanation. 
The respiratory chain generates a pmf (Ap approxi¬ 
mately -0.22 V) across the inner mitochondrial mem¬ 
brane which is approximately equivalent to a Ap' of— 
21.7 kJ-moF 1 (Nicholls and Ferguson 1997) (see. Fig. 4). 
This leads to a value of AGATP.cyt in the range between - 
58 and -65.0 kJ-mol 1 under conditions of rest, gener¬ 
ated in two steps. The first step consists in the F 0 Fi- 
ATPase complex harvesting the energy of three protons 
to build one ATP molecule in the matrix space: thus the 
mitochondrial F ( )F r ATPase complex acts as a potential 
converter of Ap' in its equivalent the inner mitochondrial 
phosphorylation potential (for more details see Ferguson 
2000). This can be calculated in terms of AGATP.mio ac¬ 
cording to Ferguson (2000) and Williamson (1979) us¬ 
ing: 

10) AG ATP ,mit(kJ _I ) ~x x Ap'tkJ- 1 ) (10) 

If x = 3 is a fixed number of transferred protons 
AGATP.mit ought to be about 3 x-22 =-66 kJ-moF 1 . 
However, because of a less than optimal efficiency, the 
resulting AGATP.mit is about -55 kJ-mol 1 (Fig. 4). The 
energy of another proton is then added to AG A TP,mit at 
the level of the ATP 4 /ADP 3 translocase (Ferguson 
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Fig. 2 Phosphorylation relationships on a semi logarithmic scale. 
For every constant pH, an iso-pH-equilibrium curve for [ADP], 
[AMP], [ATP] and [PCr] as the function of the sum of 
[ADP] + [AMP] + [P]j is calculated according to Eqs. 6, 7 and 8. 
Equilibrium constants and concentrations as in Fig. 1. It can be 
seen that free [ADP] at rest is of the order of 0.004 mmol-kg , close 
to the values calculated from 3l P-NMR-spectroscopy data, to 
attain values at fatigue of about 0.2 mmol-kg 1 at a pH of about 
6.4. From rest to fatigue free cytosolic [ADP] increases by a factor 
of 60. For abbreviations see list 

2000). The resulting driving force AG ox>ap at the inner 
mitochondrial membrane is described by: 

11) AG ox , ap (k.r 1 )~(x+l) x Ap'(kJ • mol" 1 ) (11) 

Again, because of the less than optimal efficiency the 
resulting AG ox , ap amounts to about -72.5 kJ-mol 1 or 
less (Fig. 4). 

The second step consists in a fall of AG ox>ap to the 
actual value within the cytosol. The resulting potential 
difference OAG (approximately equal to AG OXjap - 
AGATP.cyt) is the driving force which is necessary to 
sustain the flow of ATP from inside the mitochondrian 
to the cytosol. Since the flow rate of ATP is also pro¬ 
portional to the difference AG ox , ap - AGATP.cyt it also 
follows that OAG is approximately linearly related with 
the difference in the rate of OxP between rest and 
F0 2ma x- Making this assumption the effect of an in¬ 
crease of free cytosolic [ADP] itself can be described by 
the MM-kinetics, associated with a second driving force 
OAGadp given by the difference: 

12) 0AGadP= kSlDG X (AGox.ap — AGATP.cyt) (12) 


between the inner mitochondrial membrane and the 
cytosol. Using this relationship Eq. 9 can be rewritten 
under these restrictions as: 

13) F0 2 , ss =F0 2 , max /{T + k Sl /([ADP] x 0 AG A dp)} (13) 

This is consistent with the known structure of the gen¬ 
eration of mitochondrial ATP, including two indepen¬ 
dent but tightly coupled processes: firstly, the generation 
of the Ap which in turn generates AG ox , m j, by the forma¬ 
tion of ATP inside the mitochondria and, secondly, the 
delivery of ATP from inside mitochondria to the cytosol. 
There are many facts that support this assumption: 

1. Mitochondria in vitro exhibit only a slight decrease of 
the apparent Ap with an increase of the rate of V0 2 
(Nicholls and Ferguson 1997) which is assumed to be 
reflected in AG ox>ap remaining also approximately 
constant. That means that the rate of OxP increases 
linearly in relation to the difference AG ox , a p AGATP.cyt 
between the rate at rest and the maximal rate (see 
Figs. 1, 3 and 4). 

2. In resting well oxygenated muscles with a low rate of 
OxP the difference between AG ox>ap and AGATP.cyt 
comes close to 0 and the system comes close to 
equilibrium (see Figs. 1, 3 and 4). 

In Eq. 13 an equivalence between the change of 
[ADP] and the change of AGATP.cyt is proposed. Calcu¬ 
lated using Eq. 5 and 6 at pH m = 7.0 a tenfold increase 
of free [ADP] from 0.01 to 0.1 mmolT 1 would be as¬ 
sociated with a decrease of AG A TP,cyt of about 8 kJ-mol 1 
towards 0. But this would mean also that an increase of 

































325 



30 

28 

26 

24 

22 

20 

is 5 1 

O 

16 E 
_E 

14 O 
CL 

12 q: 

i- 

10 < 
8 
6 
4 

2 


0.001 0 . 


002 


0.005 

J 


0.01 0.02 


0.05 


0.1 


0.2 


0.5 


1.0 


2.0 


Rest 


ADP (MHol/kg) 


Y 

Fatigue 


Fig. 3 Steady-state activities of oxidative phosphorylation (OxP 
being approximately equivalent to FCH.ss) and glycolysis (approx¬ 
imately equivalent t oVLa), expressed as percentages of FC> 2 max and 
FXfl max , as functions of log free cytosolic [ADP J. The thick lines refer 
to values at pH = 7.0. The MM-kinetic having a 50% activity 
constant Km(0.5) of 0.027 mmol-l 1 free [ADP] is indicated by an 
arrow. The activation curve according to Eq. 13 is based on a 50% 
activity constant ksj of 0.035 mmol-l 1 free [ADP] at pH = 7.0. The 
VO 2 .SS curve which represents the rate of OxP and the curve of [PCr] 
are mirror images displaying a linear relationship within the 
physiological range of OxP activity. The cytosolic [H h ] contribute 
to the difference AGoxConst-AGATpcyt during rest and moderate 
activation of the respiratory chain. Therefore, different iso-pH 
activation curves of OxP are calculated. Free energy AG ATP cyt 
decreases linearly with increasing log free [ADP], Activation of 
glycolysis is calculated according to Eq. 16. For abbreviations see list 

OAGadp is equivalent to an increase of free [ADP] as 
calculated in Eq. 12. Therefore, [ADP] would increase 
with the 2nd power of the difference AG OXiap -AGATP,c y o 
thus explaining the experimentally observed exponent 
nox being more than 1.0. The conclusion is that the rate 
of OxP is driven not only by its substrate ADP but also 
by a potential difference of free energy along the chain of 
AG ox , ap down to AG A TP.cyt, which drives: 

1. ATP formation inside the mitochondria 

2. The ATP-ADP exchange from mitochondria to the 
cytosol and performs 

3. The down-regulation of free [ADP] at the cytosolic 
side. 

To avoid calculations of the changes of the apparent 
A p' and AG ox>ap as a function of other possible influ¬ 


encing factors without a clear mathematical basis in CS 
AG ox ap is replaced by a constant (AG OXjConst ) which re¬ 
mains uninfluenced by the rate of OxP up to F02 max . It 
is assumed that AG ox>const is close to a possible equilib¬ 
rium (AG ox>const = AG A TP,cyt) at a very low rate of OxP 
related to the proton leak. In the CS (see 
Figs. 5, 6, 7, 8, 9 and 10) AG ox , const is set to -72.5 kJ- 
mol -1 . Under these assumption the coefficient ks 1DG of 
Eq. 12 has values of about 0.65-10 4 up to -0.95-10 4 
mol -1 per 1.0 J-mol -1 if ksj is 0.035 mmol-l" 1 free ADP. 

Calculation of the steady-state activity of glycolysis 

Glycolysis is regulated mainly by the activation of PFK. 
The PFK is activated by free [ADP], [AMP] and [P] ; and 
inhibited by a high [ATP], low pH m and high concen¬ 
tration of citric acid (Krause and Wegener 1996; 
Newsholme and Start 1973). It should also be noted that 
the activity of glycolysis is regulated in part also by the 
level of phosphorylase andPK. The PFK activity de¬ 
pends on the product [ADP] x [AMP] which implies 
that [AMP] is a factor which amplifies the activation of 
PFK due to [ADP] (Krause and Wegener 1996; News- 
holme and Start 1973). The inhibition of PFK caused by 
high [ATP] is not removed with the decrease of the 
phosphorylation state of the CHEP system in the cells 
because the change of ATP is negligible. The PFK-ac- 
tivity is regulated also by fructose-2,6-bisphosphate, but 
it seems that this factor keeps PFK active and is there- 
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1 <4-cytosol-► mitochondrial -4-matrix space -► 
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Fig. 4 Diagram of a potential scheme of A TP-synthesis at the 
inside of the inner mitochondrial membrane and of the ATP 
transfer to the cytosol via ATP/ADP translocase along the 
potential differences generated from AGo X ap of about -66 kJ- 
mol 1 down to AGArpcW of about -62 kJ-moP 1 at rest. The 
AG’atp'wP is assumed to be -55 kj-mol, built up by the F 0 Fi- 
complex harvesting the energy of three protons by changing the 
conformation of the ft subunit of the F 0 F r ATPase from open (O') 
to low ( L ) to tense ( T ) state. According to Klingenberg (1980) the 
ATP 4 /ADP~ 3 translocase forms an ATP-carrier complex by 
adding about -12 kj-mol -1 to AG ATP mit powered by the trans¬ 
membrane potential Ai j/ of about -160 mV. This results in a 5:1 
preference of ATP efflux from the matrix space compared to the 
ADP uptake, contributing to the down-regulation of free cytosolic 
[ADP] during rest (Klingenberg 1980). The ATP-carrier complex is 
assumed to represents AG ox ap at the cytosolic side of the inner 
mitochondrial membrane. With an increasing rate of ATP- 
consumption by the contraction process AG AT pcyt decreases, 
whereas AG 0 x a P remains approximately constant. Free cytosolic 
[ADP] increases together with the difference AGox ,ap-AG AT pcyt = 
<j>AG (see Eq. 13). For abbreviations see list 

fore prone to be regulated, or that it acts to enhance the 
magnitude of the rate control (Krause and Wegener 
1996). Using [ADP] and [AMP] as activating modulators 
of PFK-activity, the Hill equation (Newsholme and 
Start 1973): 

14) tVss = Ida,max/{1 + ks 2 /([ADP] X [AMP])} (14) 

describes the rate of lactate (pyruvate) formation at 
“steady state” in the case of a pH m above 7.4. The time 


constant of this reaction is in the order of a few milli¬ 
seconds. Therefore every change of the rate of glycolysis 
is essentially at steady state (Mader et al. 1983; Mader 
and Heck 1986, 1996). The coefficient ks 2 = 0.00035 is 
the apparent 50% activation rate constant of glycolysis 
due to PFK activation, which is achieved for a free 
[ADP] of 0.15 mmol-1 '. It should also be pointed out 
that the V| a max expressed in terms of millimoles per 
second per kilogram is a function of the concentration of 
glycolytic enzymes. 

The maximal PFK activity is achieved at a pH m 
higher than 7.2 (Danforth 1965). The rate of glycolysis 
decreases with the fall in pH m . The in vitro inhibition of 
PFK by decreasing pH ra shows a sigmoidal shape 
(Trivedi and Danforth 1966), so that the rate of lactic 
acid production falls below 10% of its possible maxi¬ 
mum at a pH m of about 6.2 (Danforth 1965; Spriet et al. 
1987a). At a normal pH m of 6.9-7.0, glycolysis can be 
activated to about 90% of its maximum (Danforth 
1965). Many studies have reported an almost complete 
cessation of glycolysis in vitro and in vivo after maximal 
repeated contractions under anaerobic conditions at a 
pH m of about 6.4 to 6.3 associated with the cessation of 
contraction (Kent-Braun et al. 1993; Kushmerick 1987; 
Taylor et al. 1986; Thompson et al. 1995). Similar be¬ 
haviour has been observed in muscle biopsy specimens 
from human muscle taken immediately after exhausting 
exercise (Sahlin et al. 1975; Spriet et al. 1987a, b). As the 
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Fig. 5 Computer simulation of a maximal load leading to fatigue 
as a consequence of PCr-depletion and acidosis caused by lactate 
accumulation in muscles of low oxidative and medium glycolytic 
power [F0 2max = 158 ml-min '-kg 1 (about 3.75% Mit Vol, and 
approximately equal to 4.0 ml 0 2 min '-kg '), average power 
output about 16.9 Wkg '. Additional time constant tTO 2 = 9.0 s, 
FXa max = 1.0 rnmol-s -kg 1 . The transient response of V0 2 , the 
rate of glycolysis VLa ss pH ( = VLa), the concentration of lactate in 
muscle ( CLa M ) and blood ( CLa B ) are similar to measurements in 
experiments. The available VLa mlix pH is related to the intracellular 
pH. The continuous decrease of VLa ss pH ( = VLa) along with the 
accumulation of lactate is compensated by the slight increase of 
free [ADP], which activates OxP. The slight but continuous 
decrease of PCr after the initial fall is compensated in part by the 
shift of the ATP-PCr equilibrium towards [ATP], which allows a 
near complete dephosphorylation of [PCr] with less influence on 
ATP-content. The stage of fatigue is reached at the complete 
breakdown of [PCr] and the occurring decline of the ATP 
concentration. This occurs when the actual VLa ( = VLa ss pH) 
comes close to VLa max pH. For abbreviations see list 


of mixed red oxidative and white glycolytic fibres is 
assumed (Spriet et al. 1987a). 

A decreasing pH m contributes to the inactivation of 
glycolysis by the shift of the ATP-PCr equilibrium to 
ATP, which decreases [ADP] and [AMP] as was the case 
for [PCr]. The Eqs. 14, 15 or 16 do not explain the 
special details of the mechanism of the regulation of 
PFK activity. They simply provide a basis for a precise 
calculation of the experimentally observed range of the 
rate of glycolysis between conditions of rest and maxi¬ 
mal activation due to a maximal power output of the 
contraction process (Mader et al. 1983; Mader and Heck 
1986, 1996). 

Steady-state characteristic of lactate elimination 
by combustion 


inhibition of PFK by increasing [H + ] is noncompetitive, 
the decrease in the function in those enzymes which re¬ 
mains prone to activation, and hence the maximal rate 
of lactate formation can be calculated as follows (Mader 
et al. 1983; Mader and Heck 1996): 

15) r lamaxP H = i V max/(l + [H+] 3 /ks 3 ) (15) 


Combining Eqs. 14 and 15: 


16 ) ^la,sspH 

— ^la,max /{(l + [H + ] 3 /ks 3 ) 

x[l + ks 2 /([ADP] x [AMP])]} 


(16) 


where ks 3 is about 10 20 2 . 

Thus under conditions of rest the formation of 
pyruvate or lactate is equal to 0.03 mmol-kg'-min 1 or 
less if an average V] a>max ) of 1.0 mmol-s -kg 1 in muscles 


Lactate oxidation in the mitochondria of working 
muscles depends on F0 2 and pyruvate (lactate) con¬ 
centration at the level of PDH (Donovan and Brooks 
1983; Katz et al. 1989; Parolin et al. 1999). It should also 
be noted here that the reaction between lactate and 
pyruvate is always close to equilibrium. The average 
amount of pyruvate oxidized per unit 0 2 consumed is 
0.01475 mmol-ml 1 0 2 (Donovan and Brooks 1983; 
Mader and Heck 1986; Mader 1999). Therefore a 
pyruvate formation of 0.03 mmolmin '-kg 1 needs 
2.0 ml 0 2 ’min -kg _1 for its combustion. The V0 2 in the 
muscles at rest is at about this level or slightly higher. 
Experimental data show that lactate elimination by mi¬ 
tochondrial combustion i’i a , ox is a linear function of the 
current F0 2 and nonlinearly related to lactate concen¬ 
tration (Donovan and Brooks 1983; Mader and Heck 
1986); it can be calculated using the following equation: 
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Fig. 6 Computer simulation of a work load corresponding to the 
anaerobic threshold establishing a maximal oxidative steady state 
at about 63% of F0 2ma x- After a small initial rise, CLa(M) and 
CLci(B) become constant, amounting to about 4 mmol-kg 
Average mechanical power is about 9.10 W-kg , other data are 
as in Fig. 5. The contribution of VLasspH (= VLa) accounts for 
about 6% of the total energy demand. More than 95% of the 
produced pyruvate is oxidized as mitochondrial fuel. For abbre¬ 
viations see list 

17) Dia.ox = 0.021 x K0 2 , a /(l+k cLox /[la-] 2 ) (17) 

where [la ] is the lactate concentration. The factor 0.021 
results from 0.01475/0.75 where 0.75 is the average 
muscle H 2 0 space. 

The k elox is equivalent to [la] of about 
a/ 2.0 = 1.45 mrnoM -1 lactate and a metabolic RQ of 0.85 
(Mader 1999); it reflects the affinity of PDH to pyru¬ 
vate as the substrate (Mader and Heck 1986; Mader 
1999). As can be seen from Eq. 17, lactate oxidation 
depends strongly on V0 2 which determines [la] ss at 
submaximal workload (Mader and Heck 1986; Mader 
1999). During recovery from a metabolic acidosis 
caused by lactate accumulation, about one-third of the 
eliminated lactate is resynthesized to glycogen which 
can be - for reasons of simplicity - added to the rate of 
lactate oxidation. 


Feedback control of the rephosphorylating 
metabolic activity 

The phosphorylation state of the CHEP system is given 
by the sum [GP] = [ATP] + [PCr]. The dynamic change 
per unit of time can be expressed by d[GP]/dt (millimoles 
per second per kilogram) as given by the difference be¬ 


tween the GP-production rate (i.e. the instantaneous 
v GP,F02,a) and v GP Vjla ) and the rate of dephosphoryla¬ 
tion (vcp.rest + v gp,e) both expressed in terms of (milli¬ 
moles per second per kilogram). Therefore, at any time, 
the rate of change of the [GP] can be expressed by the 
following non-linear differential equation: rate of change 
of [GP] equals rate of resynthesis from respiration plus 
glycolysis, minus rate of dephosphorylation from resting 
metabolic rate plus mechanical power output, or, 

18) d[GP]/dt = f Gp p o 2 , a + ^GP.u.la — (^GP.rest + ^GP,.e)- 

(18) 

In Eq. 18, I'GP.rest is the basic rate of dephosphorylation 
which is necessary to maintain life whereas vgp,£ is set 
by the energy demand of the contraction process. The 
v GP,ko2,a is the instantaneous rate of OxP including a 
time delay resulting from adjustment of the oxygen 
transport to the consumption rate. The rate of change of 
V0 2 , a can be calculated using the following differential 
equation: 

19) dF0 2 , a /df= k F02 x ( V 0 2iSS - F0 2 , a ) (19) 

(Mader and Heck 1996) where the term V0 2 ss is cal¬ 
culated using Eq. 13 and kjz 0 2 = l/xvoi is a delay rate 
constant, which is supposed to be in the range of 4-10 s 
for the skeletal muscle, leading to an apparent total time 
constant of 15 up to 35 s in the CS of the rise of V0 2 , 
approaching a steady state under conditions of a con¬ 
stant load as observed in the whole human body 
(Barstow et al. 1993; Korecky and Brandejs-Barry 1987; 
Thompson et al. 1995). In this equation the acceleration 
of the increase of V0 2 is proportional to the difference 
between V0 2 a and the proposed V0 2 , ss according to the 
current free [ADP] calculated using Eq. 13. The term 
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Fig. 7 Computer simulation of a sinus load (1.0 Hz = 60 beats- 
min -1 ) on heart muscle having a virtual F 02 max of about 
1,000 ml-min '-kg 1 relating to about 25% V mit . Metabolic power 
output is about 17% of F0 2max . Mechanical power output is 
indicated by the lower peaks. The Via and blood lactate 
concentrations are negligible. The PCr content in conditions of 
rest when no contraction is occurring is about 14.0 mmol-kg 1 
assuming a cytosolic space of about 70% of the cell volume. Heart 
muscle in humans, at the normal frequency at rest of 65 beats-min 1 
is assumed to consume about 80 ml 0 2 -min '-kg 1 in the steady 
state, with a [PCr] = 12 mmol-kg . This is about one-sixth to one- 
eighth of the maximal rate of OxP sustainable for more than a few 
minutes. The consumption of PCr per beat is quite high, the [PCr] 
content during the contraction cycle being practically constant at 
about 12.2 mmol-kg 1 . The [PCr] [ATP] quotient is 2.05, which is 
comparable to the real situation in the working heart under the 
above-mentioned conditions. For abbreviations see list 

F0 2 , a (millilitres per second per kilogram) therefore in¬ 
cludes the delayed response of V0 2 at the beginning of 
exercise, resulting from delayed oxygen transport. The 
instantaneous E0 2a is given by: 

20 ) U GP.FO2,a( mm01 • S_1 • k § _1 ) = h V02 X ^°2,a 

x (ml • s 1 kg *) 

( 20 ) 

when the coefficient b p0 2 (millimoles per millilitre) is 
1/4.3 = 0.233, because 1.0 mmol ATP or PCr is produced 
by the consumption of 4.3 ml 0 2 (corresponding to a 
P:0 2 ratio of 5.2). The glycolytic part of the GP- 
replenishment (vgp, v la) results from: 

21) PGP,!.,la (mmol • s 1 kg A = Via X Ha,ssPH 

x (mmol • s 1 kg *). 

( 21 ) 


when b v i a is assumed to be 1.4, since the formation of 
1.0 mmol pyruvate or lactate produces from 1.35 to 
1.45 mmol ATP (or PCr). The term vi ajSS-pH is calculated 
using Eq. 16. The energy demand of the contraction 
process, in watts per kilogram muscle mass, must be 
transformed into a GP consumption rate: 

22) p G p, £ (mmol/s kg) = b Pow (mmol kg 1 • s ' W ') 

x E( W) (22) 

The estimation of b Pow depends on the efficiency of 
mechanical work performed under any given set of 
conditions for muscle contraction; as such it varies 
considerably. For the human body during cycle ergom- 
eter exercise an increase of between 11.4 and 12.5 ml- 
min -1 0 2 W _1 is commonly observed. Thus, assuming 
that 4.3 ml 0 2 are required to produce 1 mmol GP, 
this corresponds approximately to 12/(4.3 x 60) = 
0.0466 mmol GP-J 1 . The reciprocal of this quantity is 
the mechanical equivalent of GP splitting, which 
amounts therefore to 21.5 kJ-moE 1 . 

The terms (v G p,no 2 ,a + VGPvia) i n Eqs. 20 and 21 are 
the given GP-production rates, which depend on the 
phosphorylation deficit of the CHEP system expressed 
either by S[A]-[ATP] or by S[C]-[PCr] which are indexes 
of [ADP] and [AMP]. Therefore substituting Eqs. 20 
and 21 into Eq. 18 one obtains: 

23) d[GP]/df = h 1()2 x IT) 2at b ; -! a x r; a ss ])lI ^ 

— bj/02 * ^^2,a,rest — bp ow X^ 

In this equation, v [a >SSjpH is the instantaneous rate of 
glycolysis, calculated using Eq. 16 as a function of the 
current concentrations of the activating factors ([ADP] 
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Fig. 8 Computer simulation of a sinus load (3.0 Hz =180 beats- 
min -1 ) on heart muscle at a power output corresponding to 70% of 
a virtual F0 2max of 1,000 ml-min '-kg ’. At 70% of the virtual 
E0 2max the peak mechanical power is 90 W-kg 1 muscle. The 
VLa mdx is assumed to be 0.5 mmol-s -1 -kg _1 , which is equivalent to 
the one sixth of the maximal oxidative power. The time constant 
for the adjustment of oxygen delivery (t s ) to the tissue is assumed to 
be about 18 s. V0 2 at rest is assumed to be 2.5 ml-min -kg 1 
muscle; Thus 70% of the total oxidative capacity represents an 
increase of 280 times the activity at rest. From the V0 2 of 80 ml 
0 2 -min 1 -kg 1 corresponding to a heart rate of 65 beat-min 1 to 
70% V' 0 2max the change in power output is about eightfold leading 
to an increase of the heart minute volume from 5 to 30 1-min '. For 
more information and list of abbreviations see text 

and [AMP]) and taking into account also the influence 
of the pH m ([H + ]). 

The rate of change of the instantaneous oxidative GP 
production can now be expressed with the aid of the 
second differential Eq. 19, which takes into account the 
time delay of the change of F0 2 resulting from oxygen 
transport: 

24) d[GP] j, 02a /dt = b y Q2 x k^ G2 

x {E0 2 , max /[l+k Sl /([ADP] (24) 
xOAGadp)] — EOf.a}} 

and where F0 2;SS is expressed as a function of [ADP] 
and ®AGadp> as given by Eq. 13. 

The differential Eqs. 23 and 24 can be solved nu¬ 
merically as time dependent functions of FCH^rest and 
E. After each time interval, the current cytosolic equi¬ 
librium concentration of the CEIEP system has to be 
recalculated using Eqs. 6, 7 and 8, including the influ¬ 
ence of a possible change of [H + ]. Therefore Eqs. 23 and 
24 contain an additional feedback from the change of 
[la ] m , the liberation of [P] ; and the variation of PC0 2 , 


which in turn results in a change of [H + ], All these 
changes have a strong influence on the PCr-ATP equi¬ 
librium, as calculated from Eqs. 6 and 7, and on the 
activity of glycolysis, as from Eq. 16. An exact calcula¬ 
tion of the resulting change in pH m as function of the 
acid-base status might be possible but complex. How- 
ever, an approximately correct linear relationship be¬ 
tween the added base or acid equivalents and the change 
of pH m can be calculated, using the the following 
equation: 

25) pH m = 7.85 + dbuff x (0.8 x [P], - [kr] m ) 

- 0.55 x logm (PC0 2 ) (25) 

where the coefficient d buff in the CS, was given a value of 
54 mval f '-pH ra . In this way Eq. 25 combines the ef¬ 
fects on pH m of the various factors mentioned above. 
These will now be discussed separately. 

It has been shown that there is a linear relationship 
between lactate concentration and pH m within the 
physiological range of pH ra (Sahlin et al. 1975; Spriet 
et al. 1987b). The PC0 2 increases with F0 2 as a result 
of increased C0 2 production. According to the 
Hasselbalch-Henderson equation, this increase has a 
direct influence on blood as well as on internal muscle 
pH m . Despite the fact that a changing PC0 2 causes a 
change in bicarbonate towards a new PC0 2 -bicarbonate 
equilibrium and in spite of the resulting delay in the 
established steady state PC0 2 at constant power output, 
the relationship between F0 2ss and PC0 2 is roughly 
described by: 

26) PC02(mmHg) = 40.0 + 55 x F0 2ia /F0 2jmax 

x (ml • min 1 kg -1 ). 


(26) 
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At an expected normal F0 2max of 120-160 ml-min _1 kg 1 
muscle, the maximal PC0 2 is calculated to be about 80 
to 95 mmHg. The relationship between the PC0 2 and 
pH m can be calculated by: 

27) pH m = 7.736 - 0.55 x log 10 [PCO 2 (mmHg)] (27) 

for the dog skeletal muscle (Saborowski et al. 1971). In 
Eq. 25, the PC0 2 dependent influence on internal pH m 
is given by the 2nd term on the right hand side of Eq. 27, 
ApH m /Alog PCO 2 = 0.55. Therefore, Eqs. 25, 26 and 27 
are part of the simulation model, to calculate the nec¬ 
essary feedback of the change of internal pH m resulting 
from [P]j liberation, the influence of lactate accumula¬ 
tion and the change of PC0 2 brought about by the 
changes of V0 2 . 

The term v G p,E( m illiinoles P er second per kilogram) 
corresponding to power output in Eqs. 18, 22 and 23 
allows us to simulate the influence of a change of 
AtjATP.cyt towards 0 on muscle strength and power 
output. The myosin ATPase is inhibited by the increase 
of free ADP and ceases to function if the cytosolic 
quotient of [ATP]/[ADP] falls below a certain level 
(Krause and Wegener 1996). A quasi-linear relationship 
between contraction force and AGATP.cyt has also been 
reported (Jeneson et al. 1995). The reduction of power 
output as function of a decreased phosphorylation state 
of the CEIEP system can be calculated as a first ap¬ 
proximation, by: 

28) F AG = cP E x (AG AT p,cyt + 39.5)/( - 39.0) (28) 

where the term cP/ is an arbitrary conversion factor 
which is assigned a predetermined value such that the 
term F AG under conditions of rest prior to exercise as¬ 
sumes a value of 1.0. During the following period of 
exercise F AG decreases according to the decrease 
AGATP.cyt, thus leading to a decreased power output 
which is larger the heavier the exercise. This mirrors the 
phenomenon of fatigue in the CS. 

The factor cP £ can also be used to determine the 
difference in the maximal ATP turnover rate of fast or 
slow myosins which is known to be different by a factor 
of 3-5. 

The AMP deaminase reaction has not been consid¬ 
ered in this paper, because this requires another set of 
differential equations which makes it more difficult to 
demonstrate that the given solution of the high energy 
phosphate equilibrium within the regulatory circuit de¬ 
scribed by Eqs. 23 and 24 is exact and satisfies the 
conditions of an accurate regulation of glycolysis and 
OxP in the cell. 


Lactate distribution in a two-compartment model 

In humans and intact animals lactate exchange between 
the lactate producing active muscle tissue (active space) 
and a portion (fraction) of the extra cellular water space 
(passive space) can be described by a two-compartment 


model assuming that the main driving force is the con¬ 
centration difference (Freund and Zoloumian 1981; 
Mader and Heck 1996). To describe the dynamic be¬ 
haviour of lactate accumulation and removal, the dif¬ 
ferential equations must include vi ajSSjpH in the active 
compartment, the net lactate exchange between the two 
compartments and the lactate elimination by oxidation 
and glycogen resynthesis. This leads to the following 
differential equations where the lactate concentration is 
calculated for each of the compartments: 

29) d[la-] m /df=-K 1 x ([la-] m - [la] b ) 

+1.35 * (ri a)SS 

,pH ^la,ox,m ) (29) 

and 

30) d[la] b /dt=F re i x [Ki x ([la“] m - [kT] b ) - n a , 0 x,b)]- 

(30) 

The [la~] m and [la _ ] b represent lactate concentration in 
the H 2 0 space of muscle and arterial blood, respectively. 
The factor 1.35= 1/0.72 accounts for the increase of 
concentration in the muscle H 2 0 space. The F rel results 
from the quotient %F m /%Fb which is the ratio of the 
active lactate producing space to the passive lactate 
distribution space. The F re i depends on the size of the 
active (% F m ) and the passive (% F b ) space, which can be 
expressed as percentage of the total body volume 
(100%). Based on experiments the available lactate dis¬ 
tribution space, i.e. the sum of active and passive lactate 
space of the human body at rest, moderate or heavy 
exercise has been found to be about 30% to 50% of the 
total body space (Mader and Heck 1986; Searle and 
Cavalieri 1972). The two compartments are taken as 
representative (but virtual) spaces, which allow an ap¬ 
proximately correct calculation of the energy yield from 
glycolysis, based on measurable concentrations time 
curves in blood and working muscles (Mader and Heck 
1996). 

Because of the lactate proton co-transport, the lower 
buffering capacity of blood and extra cellular space 
compared to muscle tissue leads to a steeper increase of 
[H + ] in the blood which leads to a decreased speed of 
lactate exchange between muscle and blood with the 
increasing blood lactate concentration. Therefore the 
rate of lactate diffusion from muscle is nonlinearly re¬ 
lated to the blood lactate concentration. As a first 
approximation, this is described by: 

31) Ki= 0.065 x [la~] b 1,4 (31) 

The coefficient 0.065 and the exponent -1.4 are esti¬ 
mated from results in experiments on humans (Mader 
and Heck 1996) and the resulting post exercise blood 
lactate concentration curves in the simulation. The whole 
body lactate oxidation (vi a>ox ) in animals or humans can 
be subdivided into v la ox m (two-thirds) and vi a ox b (one- 
third), which represent the rates of oxidative lactate 
elimination from active and passive lactate spaces, 
respectively. 
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Application of the mathematical model 

The theoretical description of the regulation of OxP and 
glycolysis summarized by the preceding equations has 
been coded in the programming language “C”. The 
differential equations are integrated by iterative solution 
using a 5th-order Runge-Kutta-Fehlberg method 
(Engeln-Miillges and Reutter 1990). The CS of the dy¬ 
namic changes of OxP and glycolysis as function of 
power output in watts per kilogram muscle were per¬ 
formed using a normal Pentium based PC. The diagrams 
yielded by the CS were plotted using the graphic soft¬ 
ware INGRAF (version 8.00. Integrated Graphics 
Library, Sugar Land, Tex., USA). In the paragraphs 
that follow, only a few examples of the simulation of the 
dynamic pattern of the metabolic behaviour will be re¬ 
ported. Numerous other simulations of exercise dura¬ 
tions, ranging from the few seconds of a sprint to longer 
lasting exercise up to 30 min (Mader and Heck 1996; 
Mader 1998), have demonstrated that there are no re¬ 
strictions in terms of time and mechanical power to the 
application of the model. 


The metabolic response of the human body 
as a function of mechanical power output 

A prerequisite for simulating the metabolic energy sup¬ 
ply as a function of the mechanical power of the human 
body is knowledge of the ratio between active muscle 
mass and total body mass. The assumption that the 
muscle mass in normal young male athletes amounts to 
35-45% of body mass, and that about 70%-80% of the 
muscle mass can be used in maximal efforts such as 
running, yields a working muscle mass of 25%-33% of 
body mass. The lactate distribution space is not equal to 
the oxygen consuming space , the latter being closely re¬ 
lated to the active muscle mass. Therefore, the coefficient 
relating average F0 2 in the working muscle mass to 
total body mass is different from that relating lactate 
space to total body space (Mader and Heck 1986, 1996). 
The transformation coefficient from body mass to active 
muscle mass therefore varies between 100%/25%=4 
and 100%/33% = 3.0. All the CS in this paper are based 
on 1.0 kg muscle mass. 

In vitro 1.0 ml F mit can consume up to 4.5 ml 
Ovmin 1 (Hoppeler 1990; Rumsey et al. 1990; Mootha 
et al. 1997). Since the normal human skeletal muscle 
contains about 3.5% mitochondria, F0 2max can be cal¬ 
culated to be 158 ml-min 1 kg 1 muscle. If it is assumed 
that the P/O-quotient is 2.6, corresponding to a re¬ 
quirement of 11.2 ml 0 2 by 2.6 mmol ATP equivalent to 
4.3 ml 0 2 for the production of 1.0 mmol ATP, the 
above calculated F0 2max is equivalent to a production 
rate of about 0.6 mmol PCr-s'-kg 1 muscle. For a 
physical education student of 70 kg body mass, and 
30% body mass as active muscle, F0 2max is assumed to 
be 3,300 ml-min '. Together with a F0 2 of about 


400 ml-min _1 at rest this amounts to 3,700 ml-min 1 
(about 53 ml-min '-kg 1 body mass). The PCr-content is 
assumed to be about 24 mmol-kg 1 muscle. The vi a , max is 
assumed to be about 1.0 mmol-s '-kg 1 which is equiv¬ 
alent to 0.5 mmol-s^-kg 1 glucose turnover during 
maximal activation of glycolysis (Spriet et al. 1987a). 

In the CS (Fig. 5) at a supra maximal load, the first 
rapid exponential response of F0 2 , up to 70% of the 
peak F0 2 , is followed by an approximately linear in¬ 
crease until the cessation of exercise. This is consistent 
with observations in many experiments and as reported 
in some papers (e. g. Barstow et al. 1993) it is difficult to 
give a plausible explanation of this phenomenon purely 
from the point of view of the experiment. However, this 
non steady state behaviour of FOi, can be easily ex¬ 
plained if it is considered to be a consequence of the 
decline of the ATP production by vi ajSSjpH , caused by the 
inhibition of glycolysis resulting from the accumulation 
of lactate. The slight continuing decrease of the phos¬ 
phorylation state of the CHEP system leads to the de¬ 
cline in power output and the increase of free [ADP]. 
The inhibition of glycolysis shifts ATP production to the 
oxidative side. The F0 2max occurs at the stage of a 
complete breakdown of [PCr]. 

In the CS of Fig. 6 an overall oxidative steady state is 
achieved at about 60% of F0 2max and low levels of steady- 
state muscle and blood lactate (Mader and Heck 1986; 
Mader 1999). The lactate steady state results from vi a>SSjpH 
and the adjustment of lactate oxidation by increasing 
[la ] m /[la ]t>, as calculated by Eq. 17 (Mader 1999). In this 
case the oxygen deficit can be calculated by a mono ex¬ 
ponential rise of the F0 2 with only a small error. In Figs. 5 
and 6 the CS are based on real changes of F() 2 rcsl up to 
peak-F0 2 , close to the F0 2max evaluated in the ergometer 
trials and related to the stage of the muscle fatigue caused 
by the depletion of [PCr] and a decline of [ATP] (Barstow 
etal. 1993,1994;Kent-Braunetal. 1993; Sahlinetal. 1975; 
Spriet et al. 1987a; Taylor et al. 1986). 


The metabolic response of muscles with an extremely 
high oxidative capacity 

It has been suggested that F0 2 and ATP production in 
muscles with a mitochondria content of more than 25%, 
such as the flight muscles of humming birds or the 
mammalian heart, in vivo, are regulated differently when 
compared to the more usual low oxidative skeletal 
muscles (Balaban 1990; Balaban and Heineman 1989; 
Heineman and Balaban 1990; Hochachka and Matheson 
1992; Hochachka and McLelland 1997). As for the heart 
muscle, the CS in Fig. 7 shows that the PCr content at 
the steady state does not change, but oscillates slightly 
during a cycle of systole and diastole at the normal fre¬ 
quency of the human heart at rest (65 beats-min , F0 2 
about 8 ml-min 1 -lOOg 1 or 17% of the observed 
F0 2max ). This agrees with the results of ! 1 P-NMR- 
spectroscopy of the heart muscle in animal experiments 
(Kantor et al. 1986). However, the speculation that the 
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“lack of variation of high-energy phosphorylated species 
suggest a tight coupling between energy production and 
consumption, with the maximum ATP-production dur¬ 
ing peak systole” (Kantor et al. 1986) is not a crucial one. 

The CS calculated relation of [PCr]:[ATP] = 2.05 
during normal activity at rest is nearly the same as that 
observed in the heart muscle in vivo (about 1.9 ±0.1) in 
animal experiments. However, the corresponding con¬ 
tent of free [ADP] has been found to be five times lower 
in the CS compared to the experimental data (Balaban 
and Heineman 1989; Balaban et al. 1986; Kantor et al. 
1986; Katz et al. 1988a, 1989; Lewandowski et al. 1987). 
This high level of free [ADP] is mainly due to a large 
excess of free [Cr] compared to the assumed [P]; 
(Balaban et al. 1986; Balaban 1990; Katz et al. 1988b, 
1989). For example, in one of the key papers (Katz et al. 
1989) the total content of creatine compounds is deter¬ 
mined to be 27.3 mmoll -1 H 2 0, and [ATP] is assumed to 
be 6.6 mmolT '. With a [PCr]:[ATP] ratio of 2.0 under 
normal activity at rest the calculated [PCr] is 
13.65 mmoll 1 . Thus an equal amount of [Cr] is un- 
phosphorylated. This excess of free [Cr] could not be 
phosphorylated because of the assumed very low free [P]; 
(about 1.2 mmol l -1 ) leading to a [PCr]:[P]j ratio of 12.5. 
Using a free [P] ; of about 1.5 mmolT 1 ([PCr]:[P] ; of 
about 9) in the CS, then 17% of F0 2max can be achieved 
with free [ADP] of about 0.01 mmol l -1 together with the 
reported AGATP.cyt of about -61 kJ-mol 1 (see Fig. 3). So 
free cytosolic [ADP] in the heart is kept very low under 
normal functional conditions. 

The maximal rate of V () 2 of 1.0 ml F mit from heart 
muscle, as measured in a suspension of about 0.01 ml 
Unit-ml -1 H 2 0, amounts to about 4.5 ml ' min 1 Ovml 1 
F m it (Mootha et al. 1997). Extrapolation of this value to 
the intact heart on the basis of a mitochondria content 
of about 22% leads to a virtual (but not real) F0 2max of 
100 ml Ovmin I 100g 1 (or 1,000 mlmin -1 -kg -1 ). This 
constitutes a speculated 300-fold increase in the rate of 
OxP from cardiac arrest with an assumed V0 2 of 
3.0 ml-min '-kg 1 to maximal power output (Hochachka 
and Matheson 1992; Hochachka and McLelland 1997). 
However, the real changes of E0 2 are presumably con¬ 
siderably lower. 

In human hearts between normal frequencies at rest 
of 65 beats-min 1 and at maximal of 190 beats-min 1 the 
F0 2 increases about eightfold from 8.0 to about 
65 ml-min 1 -1 OOg 1 . Similarly, in the canine heart the 
myocardial V0 2 at rest was measured at about 4.1 ml 
Ovmin ^lOOg 1 (Katz et al. 1989). The observed maxi¬ 
mal increase in E0 2 was about 12-fold, i.e. up to about 
50 ml Ovmin 1 • lOOg 1 . This suggests that the actual 
F0 2max of the heart, related to the mitochondria density 
in the cardiac myocytes, is dictated by the principle of 
biological similarity, as described by: 

32) F0 2 , max (ml • min -1 • ml -1 Knit) = 4.5 * %V m ° 15 . 

(32) 


On the basis of this equation, the maximal rate of 
OxP of the myocardium which can be maintained for 
2 min is reduced from 4.5 to about 2.8 ml 0 2 -min 1 • 
ml 1 F m i t . This is about 60% to 65% of the expected 
(virtual) F0 2max and agrees well with the upper limit of 
the experimentally obtained values. The metabolic rate 
of a non-contracting heart muscle was experimentally 
determined to be approximately three times the meta¬ 
bolic rate of the soleus muscle at rest (Loiselle 1987). 
This puts the limit of the maximal possible changes of 
the metabolic rate of OxP from cardiac arrest to F0 2max 
to a maximum of 100-fold. 

The CS of Fig. 8 shows the well-known transient 
“staircase behaviour of the contraction amplitude” of 
the heart muscle which is usually observed in experi¬ 
ments on isolated hearts from a stage of low activity, or 
cardiac arrest, to that approaching a steady-state con¬ 
traction amplitude (Mast and Elzinga 1987) at about 
70% of the (virtual) F0 2max . Contraction force and 
power output are calculated as a function of AGATP.cyt 
using Eq. 28. The transient behaviour of the model 
parameters is as complex as in the real heart. The time 
constant iy 02 for a delay of OxP calculated using 
Eq. 19 or 24 is set to 1.0 s, and is due to intracellular 
adjustments of OxP. The real time constant for estab¬ 
lishing a sufficient oxygen delivery in the tissue under 
these conditions is about 20 s (Korecky and Brandejs- 
Barry 1987). This is simulated by an exponential 
increase of the F0 2max>ap after the start of exercise, 
approaching real F0 2max calculated from Eq. 32, as 
described by: 

33) F0 2 , max , ap - V 0 2 , max x (1 - e - '/ 18 ). (33) 

The steady state of V0 2 and of contraction amplitude 
in the CS is attained after more than 60 s. A transient 
decrease of [PCr] to a low level of 2.0 mmol-kg 1 is also 
observed, together with a transient high activation 
of glycolysis (about 65% of v lamax and equal to 
0.5 mmol-s-kg ) which keeps contraction force rela¬ 
tively high, and disappears with the relatively slow in¬ 
crease of E0 2 to the steady-state E0 2 . Glycolysis 
oscillates with the cardiac cycle to about ±3% of the 
average remaining activation (about 10%) at the steady- 
state OxP (see Fig. 9). In Figs. 8 and 9 the CS calculated 
ratio of [PCr]:[ATP] below 1.0 is related to about a 
tenfold increase of V0 2 compared to normal activity at 
rest in canine hearts (Loiselle 1987). A ratio of 
[PCr]:[ATP] of about 1.2 is observed in rat hearts under 
normal working conditions (Katz et al. 1988a), where 
not much of a further increase of V0 2 can be expected. 
This agrees also with the behaviour of hearts in vitro 
(Kantor et al. 1986; Katz et al. 1989; Lewandowski et al. 
1987). 

The [PCr]:[ATP] ratio under external pacing and 
adrenaline stimulation in real hearts varies from 1.7 to 
1.3, being essentially constant for each experimental 
animal (Balaban et al. 1986; Heineman and Balaban 
1990; Katz et al. 1989). In the CS a twofold increase of 


Hemmingsen (1960). 
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Fig. 9 Selective enlargement of part of Fig. 8 (60-65 s, 3.0 Hz 
repeated contraction, sinus load). In spite of an approximately 
90 W-kg 1 peak power output per beat the PCr-content remains 
approximately constant at 5.0 mmol-kg '. Glycolysis exhibits 
oscillations of about one-third of the average activity of 12% of 
VLa max , which contributes to the buffering of [PCr] during the 
contraction relaxation cycle. From normal to maximal activity the 
[PCr] decrease is 7.0 mmol-kg . For list of abbreviation see text 

V0 2 leads to a [PCr]:[ATP] ratio of 1.75 corresponding 
to an established steady-state V0 2 of about 30% of a 
virtual V0 2max , related to 26% V mit . This is within the 
range of the observations in experiments (Balaban et al. 
1986). At about 40% of the assumed virtual V0 2max the 
steady state ratio of [PCr]:[ATP] is about 1.6 which in¬ 
dicates that in a medium range of OxP the CS yields a 
high phosphorylation state. Thus, it appears that a gap 
remains between the CS and the real behaviour of hearts 
in vivo during adrenaline stimulation. This is probably 
due to the inotropic effect, which is not considered in the 
presented model. However, this is consistent with the 
basic assumptions of the model and needs only to be 
implemented. Indeed, if it is assumed that adrenaline 
leads to an increase of Ap and of its equivalent AG ox>ap , 
this effect might explain the observed difference of 
the [PCr]:[ATP] ratio between hearts driven in vivo un¬ 
der adrenaline stimulation and the CS. In the low range 
of free [ADP], below 50% of F0 2max , according to 
Eqs. 9 and 13 the effect of an increase of AG oxap of 
-5,0 kJ-mof 1 is equivalent to a five- to eightfold increase 
of free [ADP]. That this may indeed be the case is also 
suggested by experimental data, which has shown that 
the NADH/NAD + redox state can alter the rate of OxP 
of isolated mitochondria and that the increase in sub- 
maximal respiratory rate induced by pyruvate excess can 


occur without large alterations in free [ADP], [P] ; and 
[ATP] (Barstow et al. 1993; Koretsky and Balaban 
1987). 

Finally, in Fig. 10 it is demonstrated that the theo¬ 
retical scheme developed above applies also to the 
minimal conditions of the CHEP system of living skel¬ 
etal muscle cells. Using creatine analogues, the PCr 
content of muscle cells is reduced to 1/3-1/4 of the 
normal values (Korecky and Brandejs-Barry 1987; 
Meyer 1988). The behaviour of the experimentally ob¬ 
tained transients of OxP and [PCr] from rest (about 
6.0 mmol-kg ’) to steady state (about 1.8 mmol-kg ') 
and in recovery (Meyer 1989) can be re-simulated very 
precisely. 


General aspects of the regulation of OxP 

The failure to provide a satisfactory explanation of the 
regulation of OxP, as deduced from experimental re¬ 
sults, is to some extent due to a lack of consistent 
quantitative calculations of known regulating factors 
and functions. Data from 31 P-NMR spectroscopy 
show that the values of free [ADP], within the con¬ 
centration range in which OxP and glycolysis are 
proposed to be regulated, are far below the 
0.2 mmol-kg 1 wet mass which is the limit for a direct 
determination in muscle biopsy samples and in 
^'P-NMR spectroscopy (Arnold et al. 1984; Connett 
and Honig 1989; Kushmerick 1987; Kushmerick et al. 
1992; Sahlin et al. 1975). Consequently, one of the 
major sources of the discrepancies in the interpretation 
of the experimental data in the literature arises from 
the calculation of free [ADP], 
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Fig. 10 Computer simulation of animal experiments using creatine 
analogues for reducing [PCr] to 25% of the normal PCr content. 
The experiment is reported at the stage of 6 weeks using creatine 
analogues and is recalculated. With the exception of a small 
transient overshoot the simulated parameters fit the results of the 
experiment. For abbreviations see list 


Indeed, as explained above, it is difficult to make 
exact calculations of free [ADP] in this range of con¬ 
centrations from measurements of [Cr] and [ATP] and of 
the ratio of [ATP] to [PCr], as has commonly been done 
(Balaban et al. 1986; Conley et al. 1997; Connett and 
Honig 1989; Heineman and Balaban 1990; Katz et al. 
1989; Kushmerick 1987; Kushmerick et al. 1992). Free 
[Cr] is determined frequently to be in large excess, 
compared to free [P] ; . This excess of [Cr] over [P] ; adds a 
relatively large base source of calculated free [ADP] to 
the possible functional changes of free [ADP] which are 
restricted to values of [P] ; greater than 0. The calculated 
additional changes of free [ADP] are also rather small 
and seem negligible (Balaban 1990; Balaban and 
Heineman 1989; Connett and Honig 1989; Heineman 
and Balaban 1990; Katz et al. 1989). If, however, the 
calculation is based on free [P]j (Barstow et al. 1994; 
Chance et al. 1985a, b) or with only a small excess of [Cr] 
over [P] ; (Arnold et al. 1984; Kushmerick et al. 1992; 
Nioka et al. 1992; Taylor et al. 1986), then small changes 
in free [Cr] will cause large changes in the content of free 
[ADP] calculated from Eq. 2. In this case, the calculated 
change of free [ADP] in vivo fits well with the known 
Km(0.5) = 0.035 ±0.01 mmol-1 1 for the 50% activation 
of VO 2 (Arnold et al. 1984; Barstow et al. 1994; Chance 
et al. 1985b; Kushmerick et al. 1992; Nioka et al. 1992; 
Taylor et al. 1986). Thus, it seems that the only way to 


overcome this problem is to calculate the PCr-ATP 
equilibrium using free [P]j under in vivo conditions in¬ 
stead of free [Cr]. 

One of the main points of disagreement stems from 
the assumption that free cytosolic [ADP] drives OxP as 
a key substrate according to the MM-kinetic model. 
This is a difficult question because free cytosolic [ADP] 
is separated from the ATP-generating F 0 /F r ATPase 
complex by the inner mitochondrial membrane; there¬ 
fore only free [ADP] mit together with [P];, m i t inside the 
mitochondrial space are the key substrates for ATP- 
production (see Fig. 4). The [ATP] mit :[ADP] mit ratio is 
not greatly dependent on pH in the matrix space, so 
that pH does not affect this ratio, according to Eq. 4. 
According to Eq. 10 it is assumed, that the F 0 /F r 
ATPase complex generates a A(7 AT p >mit below 3 x - 
22 = -66 kJ-mol, because the reaction is potentially 
irreversible. However, the [ATP] mit :[ADP] mit ratio at 
state 4 and 3 respiring mitochondria was experimen¬ 
tally determined to be 10 up to 20 times lower com¬ 
pared to the cytosol (Erecinska and Wilson 1982; 
Klingenberg 1980; Soboll 1995). If the cytosolic and the 
mitochondrial ATP:ADP-ratio differ by a factor of 10 
or more, this rules out the possibility that OxP is di¬ 
rectly driven by free cytosolic [ADP]. The same applies 
to the exchange of ADP compared to ATP by the 
adenine nucleotide translocase which in this case is also 
not driven by concentration differences of free [ADP] 
(Erecinska and Wilson 1982; Klingenberg 1980; Soboll 
1995; Williamson 1979; Wilson et al. 1983). The only 
conclusion therefore is that a difference of potentials 
rules the rate of OxP and determines free [ADP] in the 
cytosol. If this is so, AG ox , ap which appears at the cy¬ 
tosolic side as driving force related to the Ap' must be a 

















336 


composite of a low AGATP.mit resulting from a low 
[ATP] mit :[ADP] mit ratio and an additional driving 
force generated at the level of the ATP^/ADP 3 
translocase (Ferguson 2000). As the ionic exchange by 
the ATP 4 /ADP 3 and OH'/Pi 2 translocators together 
with H + is 0 an additional -22 kJ-moF 1 is available to 
overcome the difference of a lower A(7 A i i>.mii against 
AGatp.cvi and to drive ADP transfer to mitochondria 
against the concentration difference of low free cyto¬ 
solic [ADP] to higher [ADP] mit (see Fig. 4). That the 
rate of ADP-ATP exchange across the inner mitoc¬ 
hondrial membrane can be modified or driven by 
potential differences was demonstrated by Klingenberg 
(1980). 

On the other hand it is unlikely that [ADP] mit is high 
at rest so that it does not change more than four times 
between OxP at rest (about 0.05 mmoll L ) and ITUmax 
(0.2 mmol l 1 ). Furthermore, it is also unlikely that the 
FoFj-ATPase complex is maintained always at a high 
saturation level by its substrate ADP - which indeed is 
the case - but that a change of Ap' cannot influence 
AGATP.mit- This last point is indeed at variance with the 
experimental observation that AGATP.cyt can achieve up 
to -70 kJ-moF 1 during recovery from heavy exercise 
calculated from 31 P-NMR spectroscopy experiments 
(Wackerhage et al. 1998). Since the calculations of a low 
[ATP] mit :[ADP] mit ratio are in all cases based on chem¬ 
ical determinations of [ADP] mit after rapid deprotein- 
ation and clamp freezing, the inner mitochondrial 
phosphorylation state may be deteriorated by this pro¬ 
cess. The chemically determined [ADP] mit contains a 
tightly bonded ADP/P; fraction from the F 0 /F r ATPase 
complex, amounting to 14% of the protein content of 
the inner mitochondrial membrane. The chemically ob¬ 
tained results may not, therefore, represent the real 
[ATP] mit :[ADP] mit ratio. In addition, new experimental 
results indicate that the 3:1 coupling of proton to ATP 
formation is not as fixed as commonly believed (Fer¬ 
guson 2000). So the free [ADP] mit together with free 
[P]i,mit can be kept low and AGATP.mit may approach 
AGATP.cyt under conditions of rest. If, according to El¬ 
ston et al. (1998), four protons are transferred to pro¬ 
duce one ATP, Eq. 11 would indicate that AGATP.mit 
could be in the range of-60 to -70 kJ-mol _1 , i.e. equal to 
AGATP.cyt at rest. Since the Ap' is kept relatively constant 
even with an increasing rate of OxP up to V0 2max , 
AGATP.mit is also kept high. If, however, the [ATP] mit : 
[ADP] mit ratio as a function of Ap' is not so much below 
the cytosolic [ATP]:[ADP] ratio even under conditions 
of rest it is not necessary to give up the idea that OxP is 
controlled via the MM kinetic, regulated by the level of 
the free [ADP]. In this case indeed, ADP does not need 
to be transported against a high concentration differ¬ 
ence. This also was first concluded by Williamson 
(1979). 

If, however, free [ADP] mit together with [P| i mit is kept 
at a low level by a high affinity of the F 0 /F r ATPase 
complex for both substrates which mirrors Ap' at high 
AGATP.mit? then the assumption that the free cytosolic 


[ADP] directly drives the exchange rate in part at the 
level of the ATP^/ADP 3 translocase is very likely. 
Under these assumptions the regulation concept of OxP 
from the cytosolic side is simple and clear. The equiva¬ 
lent of the apparent Ap' (AG OXiap ) reflects the driving 
force of OxP from the mitochondrial side according to 
Eq. 11 and Eq. 13 is to some extent a causal but sim¬ 
plified mathematical description of this process which 
determines the degree of the activity of OxP from the 
cytosolic side. 

In this paper the model is limited to the cytosolic 
side of OxP with the assumption of a constant 
AG ox-const of -72.5 kJ-mol 1 . This is the equilibrium 
point at which the rate of OxP, as calculated by Eq. 13, 
equals 0. This first approach neglects existing changes 
of the apparent Ap' observed in experiments. At every 
stage of activity of OxP, even at rest, a deviation from 
the equilibrium related to a OAGadp greater than 0 is 
necessary to keep the system working. Indeed, a com¬ 
plete equilibrium cannot be reached because of the 
existence of a basic oxidation rate of 0.5 ml 02'inin 1 • 
kg 1 or more in normal skeletal muscles, related to the 
proton leak of the inner mitochondrial membrane and 
the heat production necessary to maintain body tem¬ 
perature at 37°C (Loiselle 1987; Rolfe and Brown 
1997). At complete equilibrium no fluxes can exist. It is 
thus clear that increasing fluxes need increasing driving 
forces and increasing substrate or activator concentra¬ 
tions which must deviate from the equilibrium status 
close to resting conditions. 

In summary the proposed concept of regulation is a 
basic approach at the simplest possible level; it can be 
taken as an attempt to integrate the two main dis¬ 
agreeing concepts of the regulation of OxP in vivo 
from the past. Results from the CS demonstrate a 
high level of similarity with the experimentally ob¬ 
served dynamic behaviour of measurable parameters 
of OxP, glycolysis and the behaviour of the parame¬ 
ters of the CHEP system. The advantage of a com¬ 
plete simulation of the behaviour of OxP and 
glycolysis as function of the energy demand of muscle 
cells is that the model also calculates from a clear 
theoretical base the details of the parameters that can 
be measured in various experimental situations. Fur¬ 
thermore, the CS shows that complex dynamic be¬ 
haviour in real biological systems can be ruled by a 
few general principles and a relatively simple set of 
equations. It is, however, clear that the model needs a 
more detailed implementation of the characteristics of 
the reactions occurring at the site of the mitochondria. 
This however does not change the calculation of the 
main regulation of OxP and glycolysis as explained in 
this paper. 
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